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ABSTRACT 

The inspiral and coalescence of double neutron star (NS-NS) and neutron star-black hole (NS-BH) 
binaries are likely to be detected by advanced networks of ground-based gravitational wave (GW) in- 
terferometers. Maximizing the science returns from such a discovery will require the identification and 
localization of an electromagnetic (EM) counterpart. Here we critically evaluate and compare several 
possible counterparts, including short-duration gamma-ray bursts (SGRBs), "orphan" optical and 
radio afterglows, and ~ day-long optical transients powered by the radioactive decay of heavy nuclei 
synthesized in the merger ejecta ( "kilonovae" ) . We assess the promise of each potential counterpart in 
terms of four "Cardinal Virtues" : detectability, high fraction, identifiability, and positional accuracy. 
For viewing angles within the half-opening angle of the jet (0 o bs ^ &j) the SGRB and associated af- 
terglow are easily detectable within the range of Advanced LIGO/Virgo if the jet energy (Ej) and the 
circumburst density (n) are similar to those inferred from existing SGRB observations. For modest 

7 /S 

off-axis angles (# b s J$ 20j), the orphan optical afterglow is detectable with LSST if £^,50 n > 0.002; 
the fraction of such events is ~ 79j ~ 0.1. At even larger viewing angles (i.e., the majority of observers) 
the isotropic kilonova emission dominates, with a peak optical brightness of ~ 19 — 22 mag within 
the Advanced LIGO/Virgo volume, detectable with LSST using a specialized 1-day cadence. Radio 
afterglow emission from an initially off-axis jet or from sub-relativistic ejecta is also isotropic, but 

peaks on a timescale of months-years; this signal is detectable provided that £^,50 n^ 8 (u/c) 11 / 4 > 0.2 
(for off-axis afterglows, v/c ~ 1). However, existing SGRB afterglows do not satisfy this criterion, 
indicating a low probability of radio detections. Taking into account the search strategy for typical 
error regions of tens of square degrees, our primary conclusion is that SGRBs are the most useful 
EM counterparts to confirm the cosmic origin of a few GW events, and to test the association with 
NS-NS/NS-BH mergers. However, for the more ambitious goal of localizing and obtaining redshifts 
for a large sample of GW events, kilonovae are instead preferred. Off-axis optical afterglows will be 
detectable for at most ~ 10% of all events, while radio afterglows are promising only for the unique 
combination of energetic relativistic ejecta in a high density medium, and even then will require hun- 
dreds of hours of EVLA time per event spread over months-years. Our main recommendations from 
this analysis are: (i) an all-sky 7-ray satellite is essential for temporal coincidence detections, and for 
GW searches of 7-ray triggered events; (ii) LSST should adopt a 1-day cadence follow-up strategy, 
ideally with ~ 0.5 hr per pointing to cover GW error regions (the standard 4-day cadence and depth 
will severely limit the probability of a unique identification); and (iii) radio searches should only focus 
on the relativistic case, which requires observations for a few months. 

Subject headings: gamma rays: bursts-gravitational waves-binaries-stars: neutron 



1. INTRODUCTION 



The first direct detection of gravitational waves 
(GWs) is anticipated within the decade once the 
groun d -based interferomete rs LIGO 4 (jAb ramovi ci et al 



1992; I Abbott et all 12009ft and Virgo 5 (|Caron et al 



1999; Acernese et al. 2009) are upgraded to "ad 



vanced" sensitivity (hereafter ALIGO/Virgo). The Large 
Scale Cryogenic Gravitational Wave Telescope (LCGT; 
iKuroda fc LCGT Collaborationll2010l ) is under construc- 
tion in Japan and is anticipated to join ALIGO/Virgo 
by about 2018. The most promising astrophysical 
GW sources in the frequency range of these detec- 
tors are the inspiral and coalescence of compact ob- 
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ject binaries with neutron star (NS) and/or black hole 
(BH) constituents. Although this accomplishment will 
stand on its own merits, optimizing the science re- 
turns from a GW detection will require the identification 
and study of coinc i dent e lectromagnetic (EM) c ounter- 
parts /e.g.. iSchutzl [19861 [20021: iSvlvestrd [2lm iStubbsl 
120081 iPhinnevi 120091 : IStamatikos et all l2009( l This is 
important for several reasons, including lifting degen- 
eracies associated wit h the inferred binary parameters 
(Hugh es~&: Hold |2003| ) : reducing the signal-to-noise ra- 
tio for a confident GW detection ( Kochan ek fc Piranl 
fl99l iDalal et all [20061 : lHarrv fc Fairhurstl 1201 lh : and 
identifying the merger redshift, thereby setting the 
energy scale and allowing an independent measure- 



rameters (e.g. 


Krolak & Schutzl 119871: iChernoff & Finnl 


19981 iHolz & ] 


Hughes! 20051: iDeffavet & Menoul 120071: 


Nissanke et al.l 


2010). The potential wealth of comple- 



mentary information encoded in the EM signal is likewise 
essential to fully unraveling the astrophysical context 
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of th e event (|Phinnevl 120091: iMandel fe O'Shaughnessvl 
l2010f > . for examp le an association with specific stellar 
populations fe.g.. lFong et al.lfeOlOl) . 

Motivated by the importance of EM detections, in this 
paper we address the critical question: What is the most 
promising EM counterpart of a compact object binary 
merger? The answer of course depends on the definition 
of "most promising" . In our view, a promising coun- 
terpart should exhibit four Cardinal Virtues, namely it 
should: 

1. Be detectable with present or upcoming telescope 
facilities, provided a reasonable allocation of re- 
sources. 

2. Accompany a high fraction of GW events. 

3. Be unambiguously identifiable (a "smoking gun"), 
such that it can be distinguished from other astro- 
physical transients. 

4. Allow for a determination of ~ arcsecond sky posi- 
tions. 

Virtue #1 is necessary to ensure that effective EM 
searches indeed take place for a substantial number of 
GW triggers. Virtue #2 is important because a large 
number of events may be necessary to build up statis- 
tical samples, particularly if GW detections are rare; in 
this context, ALIGO/ Virgo is predicted to detect NS- 
NS mergers at a rate ranging from ~ 0.4 to ~ 400 yr" 1 , 
wit h a "best-bet" rate of ~ 40 yr" 1 (jAbadie et al.ll2010bt 
cf. iKopparapu et al J 12008( 1. while the best-bet rate for 
detection of NS-BH mergers is ~ 10 yr -1 . Virtue #3 is 
necessary to make the association with high confidence 
and hence to avoid contamination from more common 
transient sources (e.g., supernovae). Finally, Virtue #4 
is essential to identifying the host galaxy and hence the 
redshift, as well as other relevant properties (e.g., asso- 
ciation with specific stellar populations). 

It is important to distinguish two general strategies 
for connecting EM and GW events. One approach is to 
search for a GW signal following an EM trigg er, either in 
real time or at a post-proc essing stage (e.g.. iFinn et al.1 
119991: iMohantv et al.ll2004D . This is particularly promis- 
ing for counterparts predicted to occur in temporal co- 
incidence with the GW chirp, such as short-duration 
gamma-ray bursts (SGRBs). Unfortunately, most other 
promising counterparts (none of which have yet been in- 
dependently identified) occur hours to months after co- 
alescence 6 . Thus, the predicted arrival time of the GW 
signal will remain uncertain, in which case the additional 
sensitivity gained from this information is significantly 
reduced. For instance, if the time of merger is known 
only to within an uncertainty of ~ hours(weeks), as we 
will show is the case for optical(radio) counterparts, then 
the number of trial GW templates that must be searched 
is larger by a factor ~ 10 4 — 10 6 than if the merger time 
is known to within seconds, as in the case of SGRBs. 

6 Predicted EM counterparts that may instead precede the 
GW sig nal include emission powered by the magnetosphere of the 
NS (e.g. IHansen fe Lvutikovll200U iMc Williams fe LevjnllMllll. or 
crack ing of the NS crust due to tidal interactions (e.g. I'lYoiaet al.l 
120101 1. during the final inspiral. However, given the current uncer- 
tainties in these models, we do not discuss them further. 




Fig. 1. — Summary of potential electromagnetic counterparts 
of NS-NS/NS-BH mergers discussed in this paper, as a function 
of the observer angle, o t, s . Following the merger a ccntrifugally 
supported disk (blue) remains around the central compact object 
(usually a BH). Rapid accretion lasting < 1 s powers a collimated 
rclativistic jet, which produces a short-duration gamma-ray burst 
(J2}. Due to relativistic beaming, the gamma-ray emission is re- 
stricted to observers with 6> t, s < 9j , the half-opening angle of the 
jet. Non-thermal afterglow emission results from the interaction of 
the jet with the surrounding circumburst medium (red). Optical af- 
terglow emission is observable on timescales u p to ~ days— weeks by 
observers with viewing angles of o t, a < 29 j ( i]3.1|l . Radio afterglow 
emission is observable from all viewing angles (isotropic) once the 
jet decelerates to mildly relativistic speeds on a timescale of weeks- 
months, and can al so b e produced on timescales of years from sub- 
rclativistic ejecta ( £|3,2H . Short-lived isotropic optical emission last- 
ing ~few days (kilonova; yellow) can also accompany the merger, 
powered by the radioactive decay of heavy elements synthesized in 
the ejecta (Jp. 

A second approach, which is the primary focus of 
this paper, is EM follow-up of GW triggers. A poten- 
tial advantage in this case is that counterpart searches 
are restricted to the nearby universe, as determined by 
the ALIGO/Virgo sensitivity range (redshift z < 0.05 — 
0.1). On the other hand, a significant challenge are the 
large error regions, which are estimated to be tens of 
square degrees e ven for optimistic conf i gurations of GW 
detectors (e.g., |GJ i rsel fc Tintol HgM iFairhurstl 120091 : 
I Wen fc ChefjfeoiOt iNissanke et al J 1201 1( 1. Although it 
has been argued that this difficulty may be alleviated 
if the search is restrict ed to galaxies within 200 Mpc 
(|Nuttall fc Suttonl 12010( 1 . we stress that the number of 
galaxies with L > 0.1 L * (typical of SGRB host galax- 
ies; lBerger|[2009ll2011bH within an expected GW error 
region is ~ 400, large enough to negate this advantage 
for most search strategies. In principle the number of 
candidate galaxies could be reduced if the distance can 
be constrained from the GW signal; however, distance 
estimates for individual events are rather uncertain, es- 
pecially at that low SNRs tha t will characterize most de- 
tections ((Nissanke et al.ll2010D . Moreover, current galaxy 
catal ogs are incomplete within th e ALIGO/Virgo volume 
(e.g. iKulkarni fc Kasliwalll2009( l. especially at lower lu- 
minosities. Finally, so me mergers may also occur outsid e 
of their host galaxies (jBergcr 2010a; Kcll ev et al.l i2010). 

At the present there are no optical or radio facilities 
that can provide all-sky coverage at a cadence and depth 
matched to the expected light curves of EM counter- 
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parts. As we show in this paper, even the Large Syn- 
optic Survey Telescope (LSST), with a planned all-sky 
cadence of 4 d and a depth of rw 24.7 mag, is unlikely 
to effectively capture the range of expected EM coun- 
terparts. Thus, targeted follow-up of GW error regions 
is required, whether the aim is to detect optical or radio 
counterparts. Even with this approach, the follow-up ob- 
servations will still require large field-of-view telescopes 
to cover tens of square degrees; targeted observations 
of galaxies are unlikely to substantially reduce the large 
amount of time to scan the full error region. 

Our investigation of EM counterparts is organized as 
follows. We begin by comparing various types of EM 
counterparts, each illustrated by the schematic diagram 
in Figure [TJ The first is an SGRB, powered by accretion 
following the merger Even if no SGRB is produced 
or detected, the merger may still be accompanied by rel- 
ativistic ejecta, which will power non-thermal afterglow 
emission as it interacts with the surrounding medium. In 
S|3]we explore the properties of such "orphan afterglows" 
from bursts with jets nearly aligned towards Earth (op- 
tical afterglows; £|3.1[) and for larger viewing angles (late 
radio afterglows; §3.2j) . We constrain our models using 
the existing observations of SGRB afterglows, coupled 
with off-axis afterglow models. We also provide a realistic 
assessment of the required observing time and achievable 
depths in the optical and radio bands. In £|4]we consider 
isotropic optical transients powered by the radioactive 
decay of heavy elements synthesized in the ejecta ( "kilo- 
novae" ) . In Sjl] we compare and contrast the potential 
counterparts in the context of our four Cardinal Virtues. 
Although some of these counterparts have been discussed 
previously in the literature, we examine them together to 
better highlight their relative strengths and weaknesses. 
Drawing on the properties of the various counterparts, 
in fj6] we make specific recommendations for optimizing 
the follow-up with 7-ray satellites, wide-field optical tele- 
scopes (PTF, Pan-STARRS, LSST), and radio telescopes 
(EVLA, ASKAP). We summarize our conclusions in SJ71 

2. SHORT-DURATION GAMMA-RAY BURSTS 

The most commonly discussed EM counterpart of NS- 
NS/NS-BH mergers is an SGRB, powe red by accretion 
onto the central co mpact object (e.g., [P aczvnsk il 119861: 
Eichler et all Il989t iNaravan et "all ll992nRezzolla et all 
201 1|) . The Swift satellite, and rapid follow-up observa- 
tions with ground-based telescopes, have revolutionized 
our understanding of SGRBs by detecting and localizing 
a sign ificant number of their afterglows for the first time 
(e g IBerger et all 120051: iFox et al.l 12001 iHiorth et all 
120051 Bloom et all \200&) . This has enabled the dis- 
covery that SGRBs originate from more evolved stel- 
lar populations than those of long-duration GRBs, con- 
sistent with an origin associated with NS-NS mergers 
(IBerger et al 1120051 iBloom et al l[2006l [Leibler fc Bergerl 
I20ld iBergerl l2011ri iFong et al.l l2011al) . The study of 
SGRB afterglows has also established a scale for the en- 
ergy release and circumburst density that are lower than 
for long GRBs, with E < 10 51 erg and n < 0.1 cm" 3 
(|Berger et al.ll2005l:ISoderberg et al.ll200llBergerl2007al ). 
These observations have also provided evidence for col- 
limation in at least one case (GRB 051221A), with a jet 
half-opening angle of 9j ~ 0.12 (|Burrows et al.l 120061 : 
ISoderberg et al.ll2006D . and upper or lower limits in ad- 
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Fig. 2. — Cumulative detection rate of SGRBs with measured 
redshifts > z (thick solid lin e), calculated using 19 (mostly Swift) 
SGRBs (e.g., Bcrgcr 2011b). Dashed vertical lines mark the esti- 
mated sensitivity range of ALIGO/Virgo to NS-NS and NS-BH 
mergers, respectively, including a boost due to the face-on bi- 
nary orientation. The thin solid line shows an approximate fit 
to Nqrb obs(> z ) at low redshift. The dot-dashed line shows an 
estimate of the total SGRB detection rate (with or without redshift 
information) by an all-sky 7-ray telescope with a sensitivity similar 
to Fermi/GBM. 

ditional cases (|Fox et al.ll2005l : lGrupe et al.ll2006l : IBerger! 
l2007bf) . overall suggestive of wider opening angles than 
for long GRBs. 

Despite this progress, it is not yet established that 
all SGRBs ar e uniquely associated with NS-NS/NS-BH 
mergers fe.g.. iHurlev et al.ll2005l: iMetzger et~altl2008bf) , 
nor that all mergers lead to an energetic GRB. The 
energy of the GRB jet, for instance, may depend sen- 
sitively on the mass of the remnant accretion disk, 
which from numerical simulations appears to vary by or- 
ders of magnitude (~ 10~ 3 — 0.1 M©), depending on 
the propertie s of the binary and the high-density equa- 
tion o f state dRuffert et al.lll997t Ijanka et al J 119991: Eei 
20011: Rosswo getalJ 120 03: Shibata & T aniguchil I2008F 
Duez et al.ll2010HChawla et al.|[2010h . 

Although SGRBs are bright, they occur relatively 
rarely within the range of ALIGO/Virgo. To illustrate 
this point, in Figure [5] we plot the cumulative rate 
at which SGRBs are currently detected above a red- 
shift z, VGRB,obs(> z). This distribution includes 19 
SGRBs with well-determined redshifts, obt ained from 
host galaxy associations (e.g., IBergerl [20091) . Since its 
launch in late 2004 Swift has detected SGRBs at a 
rate of ~ 10 yr _1 , of which ~ 1/3 have measured red- 
shifts. Shown for comparison are the sensitivity ranges 
D r « 1.5 x 196[410] w 295[615] Mpc for detection of 
NS-NS [NS-BH] mergers by ALIGO/Virgo 7 , where the 
factor of « 1.5 (included only in this section and i j3.1[) 
accounts for the stronger GW signal from face-on merg- 
ers, which characterize the geometry of GRB jets (e.g., 

7 Throughou t this paper w e adopt the fiducial values for D r RJ 
200 Mpc from [Abadic et al. (2010b), who define detections as 
events with signal-to-noise ratio (SNR) of 8 in a single detector, 
assuming NS/BH masses of 1.4/IOMq. This choice is conservative 
because for a network of N detectors, the sensitivity range at fixed 
SNR increases D r <x A" 1 / 2 . On the other hand, the real detection 
range of a network depends on the data quality (e.g., Gaussianity 
and stationarity) and detection pipeline. Once a value for D r is 
chosen, all of the results presented in this paper may be rescaled 
accordingly. 
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iKochanek fc Piranlll99l lSchutdl20ll . 

Figure [2] illustrates the striking fact that no SGRBs 
with known redshifts have yet occurred within the 
ALIGO/Virgo range for NS-NS mergers, while only two 
SGRBs (061201 and 080905) have occurred within the 
NS-BH range. Though selection effects and low-number 
statistics undoubtedly distort the true redshift distribu- 
tion from that shown in Figure [2 at low redshift the 
distribution should nevertheless scale as -/Vgrb.oIjs °c z 3 . 
By fitting the lowest redshift bins to a distribution of this 
form, we find that < 0.03(0.3) SGRBs per year are cur- 
rently being localized by Swift within the ALIGO /Virgo 
range for NS-NS(NS-BH) mergers 8 . Thus, even assum- 
ing that Swift (or a mission with similar capabilities) 
operates simultaneously with ALIGO/Virgo, SGRBs are 
clearly not ideal counterparts to localize a large num- 
ber of mergers. Obtaining a single GW redshift in this 
fashion could require a decade of observations. 

Localization is of course only one desirable virtue of 
an EM counterpart. Due to the short duration of both 
SGRBs and the GW signal, and the short expected de- 
lay (< seconds) between them, a time coincidence be- 
tween these events is sufficient to enable a statistically 
confident association. Even if the redshift cannot be ob- 
tained, a coincident detection will still confirm the astro- 
physical nature of the GW signal, prove the connection 
between SGRBs and NS-NS/NS-BH mergers, and allow 
studies of the dependence of th e binary inclination on th e 
properties of the GRB jet (e.g. . IKochanek fc Piranl 1993). 
Coincidence searches for GW bursts using the time and 
sky coordinates of detected SGRBs were already con- 
ducted durin g previous LIGO /Virgo Sc ience Runs (e.g., 
lAbadie et alll2010at lAbbott et al.l)2010D 

To estimate how long ALIGO/Virgo must operate be- 
fore a connection between SGRBs and NS-NS/NS-BH 
mergers can be tested, we also plot in Figure [2] an esti- 
mate of the low-redshift distribution, but including all 
detectable SGRBs (with or without redshift informa- 
tion), which we estimate by multiplying the "with red- 
shift" distribution by a factor « 10. This factor accounts 
for the higher rate, ~ 20 yr _1 , that Fermi's Gamma- 
Ray Burst Monitor (GBM) detects SGRBs, as compared 
to the rate with redshift from Swift (correcting also for 
the GBM field of view, which covers only ~ 60 per cent 
of the sky). This estimate illustrates that a few Fermi 
bursts over the past few years probably occurred within 
the ALIGO/Virgo volume. Thus, an all-sky 7-ray mon- 
itor with a sensitivity similar to Fermi/GBM could test 
whether SGRBs originate from NS-NS/NS-BH mergers 
within just a few years after ALIGO/Virgo reaches full 
sensitivity, even if it does not lead to a significant im- 
provement in the sky localizations. 

One issue raised by the above analysis is that the 
observed SGRB rate within the ALIGO/Virgo volume, 
even when corrected for partial sky coverage, is much 
lower than t h e best -bet NS-NS merger rate of ~ 40 yr" 1 . 
iNakar et al.l ((2006) estimate that the local volumetric 
SGRB rate is > 10 Gpc -3 yr _1 , which corresponds to 

8 The sensitivity range for a GW detection may be increased 
somewhat if the search is restricted to the time interval and sky 
position of the SGRB in the case of a 7-ray triggered search 
IIKocha nck & Piran 1993]), but this does not alter our conclusion 
that SGRBs are a rare occurrence in the range of ALIGO/Virgo. 



an all-sky rate of Nqhb all-sky ~ 0-3 yr 1 at a distanc e 
of < £> r ,NS-NS ~ 200 Mpc (cf., iGuetta fc Piranl 12005]) . 
consistent with our estimates in Figure [2] and still two or- 
ders of magnitude below ~ 40 yr" 1 . Reconciling this re- 
maining discrepancy requires either that the true merger 
rate is lower than the best-bet rate; that all mergers are 
not accompanied by a bright SGRB; or that the 7-ray 
emission is beam ed (e.g.. lRosswog fc Ramirez-RuizfeOOa 
lAlov et al.ll2005l) . 

Expanding on this final possibility, if the typical SGRB 
jet has a half-opening angle 0j < tt/2, then only a frac- 
tion fb n w 1 — cos 9j w 1 of viewers with observ- 
ing angles o bs ^ 0j will detect a bright SGRB. For all 
other observers (the majority of cases) the prompt emis- 
sion is much dimmer due to relativistic beaming. Recon- 
ciling the "observed" and best-bet rate by beaming alone 
thus requires fb n ~ 0.01, or 0j ~ 0.12, similar to the 
open i ng angle inferred for G RB051221A (Bu rrows et al.1 
l2006HSoderberg et al.ll2006l ). 

A mystery associated with SGRBs is that ~ 1/4 — 
1/2 are followed by variable X-ray emission with 
a fluence compa rable to, or in exces s of the ini- 
tial burst (e.g., iNorris fc Bonnelll 120061: iPerlev et all 
2009). Although the origin of this extended emis- 
sion is still debated, one explanation is that it re- 
sults from ongoing energy output from a highly mag- 
netized neutron star, which survives the NS-N S merger 
(|Metzger et alJl2008bt iBucciantini et al.ll2011l ). Regard- 
less of its origin, if some mergers are indeed accom- 
panied by extended X-ray emission, this provides an 
additional potential EM counterpart, especially if the 
X-ray emission is more isotropic th an the SGRB it- 
self (as predicted by several models: | MacFadven et al.l 
20051 iMetzger et all l2008bt iBarkov fc Pozanenkol 120111: 
Bucciantini et aLll2011l ) . Considering alternative prompt 
counterparts is germane because the lifetime of Swift and 
Fermi are uncertain, while the next generation of pro- 
posed high-energy transient satellites (e.g., Janus, Lob- 
ster) are most sensitive at soft X-ray (rather than 7— ray) 
energies, which could reduce their sensitivity to detect- 
ing the prompt SGRB phase. The difficulty of relying 
on this extended X-ray signal is twofold: (i) separating 
such cases from soft long GRBs; and (ii) these events 
represent only a fraction of all SGRBs. 

3. AFTERGLOWS 

Even in the absence of an SGRB, an orphan afterglow 
may prov ide a bright electromagnetic link to a GW trig- 
ger (e.g., iCoward et aJll20Tl iNafair fc Piradl20ll . The 
orphan afterglow can be on-axis if the 7-ray emission was 
missed due to incomplete sky coverage by 7-ray satellites, 
or it can be off-axis if the relativistic jet was initially 
pointed away from our line of sight. For off-axis observers 
the afterglow emission peaks at a later time and at a 
lower brightness level than for on-axis observers, making 
the detection of a counterpart more challenging. How- 
ever, a higher fraction of events, oc 0^ hs , occur at larger 
angles, with the total fraction of detectable counterparts 
depending on the largest viewing angle at which emis- 
sion is still detectable. On a timescale of ~ days after the 
merger, the afterglow emission is still partially beamed 
and peaks at optical wavelengths ( §3.1|) . At later times, 
weeks-months, the emission is mostly isotropic and peaks 
at radio wavelengths, once the jet decelerated to mildly 
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relativistic velocities, fi < 1 (2L2J). 

3.1. Optical Afterglow 

To gain insight into the afterglow emission that may 
accompany a merger event (on- or off-axis), we use ex- 
isting observations of SGRB optical afterglows discov- 
ered in rapid follow-up observations. Figure |3] shows 
optical d etection s and u pper limits from the c ompi- 
lation of iBergerl (|2010a| ) and iFong et all (|2011al ). ex- 
pressed in luminosity and in apparent magnitude for 
a source at a distance of 200 Mpc. Shown for 
comparison are on-axis SGRB afterglow model s with 
6» obs ~ 0j from Ivan Eerten k MacFadvenl (|2011h 9 (see 
also Ivan Eerten et aJj l2010f ) that span the range of de- 
tected afterglows. Also shown are the sensitivity limits of 
existing and planned wide-field survey telescopes (PTF, 
Pan-STARRS, LSST), taking into account that about 10 
pointings will be required to cover a typical GW error re- 
gion, leading to at most 0.5 hr per pointing to cover the 
region in a single night; for LSST we also show the ex- 
pected depth of the normal survey mode (r « 24.7 mag), 
which can cover a GW error region in only a few minutes. 
The expected maximum depth is about 22 mag for PTF, 
23.5 mag for Pan-STARRS, and 26.5 mag for LSST. We 
note that the trade-off between limiting magnitude and 
localization area (A) is simply Am » — 2.51og(A 1 / 2 ) such 
that for a best-case scenario of A ~ few deg 2 , these tele- 
scopes can achieve a greater depth by about 1.2 mag in 
a single night. 

Figure [3] demonstrates that for a typical jet half- 
opening angle of dj — 0.2, afterglow models with jet en- 
ergies of Ej w 10 — 10 50 erg and circumburst densities 
of n 10~ 3 — 1 cm~ 3 are consistent with the range of 
observed optical luminosities. We can define a figure of 
merit for the combi nation of energy and density (e.g., 
iGranot k Saril l2002h . which based on the observed on- 
axis optical afterglow luminosities has an upper bound 
of: 

FOM opt , oa = nl /2 < 0.1 (1) 

and a mean value for the detected sample of 
FOM opt: on ~ 0.01; here Ej y50 is the jet energy in units of 
10 50 erg, riQ is the circumburst density in units of cm -3 , 
and we assume a typical value of p — 2.5 for the electron 
power law distribution. We note that the sample in Fig- 
ure [3] represents all SGRBs with deep optical searches, 
and hence also a detected X-ray afterglow. Since ~ 1/4 
of SGRBs lack detected X-ray afterglows, and not all 
events with X-ray detections had deep follow-up optical 
searches, it is possible that some SGRB optical after- 
glows are dimmer than those in Figure |3l leading to an 
even lower mean value of F0M op t jOn than inferred above. 
Nevertheless, we conclude that ~ 1/2 of SGRBs within 
the range of ALIGO/ Virgo (even those missed due to 
incomplete 7-ray sky coverage) should produce optical 
emission detectable by LSST for at least ~ 10 d; the 
brightest events should be detectable for a few days even 
by less sensitive surveys such as PTF. 

Using the range of energies and circumburst densities 
inferred for on-axis afterglows, we can now predict the 
appearance of off-axis light curves. Figure [4] shows the 

9 http://cosmo.nyu.edu/afterglowlibrary/ 



same range of models from Figure |3l but now for an ob- 
server angle of 6> bs = 29 j . For the parameters of existing 
SGRB afterglows, the peak flux at even larger viewing 
angles (e.g., 9 a b s = A9j; Figure [5]) is too low to be de- 
tected even with LSST. For the range of applicable mod- 
els, the off-axis light curves rise to maximum brightness 
on a timescale of ~ 1 — 20 d, with a peak luminosity 
of ~ 10 38 — 10 41 erg s _1 ; this corresponds to a appar- 
ent brightness of > 23 mag at 200 Mpc. The models 
with n = 10 -3 cm -3 peak on a timescale about 7 times 
longer, and with a luminosity that is about 300 times 
lower, than those with n = 1 cm' 3 . We also note that 
the off-axis afterglow light curves in the higher density 
cases are qualitatively similar to those of kilonovae, al- 
though the latter fade more rapidly after the peak and 
have a distinct color evolution (jJIJ. 

For the off-axis light curves we can define a separat e 
figure of merit (c.f., Equation 11 of lNakar k Pirar3l2011l) : 

FOM optiOS = E jt50 n r /8 . (2) 

We define detectable cases as those rising by at least an 
order of magnitude (2.5 mag) above the LSST maximal 
depth, corresponding to FOM opt , ff > 0.002 (Figures 0] 
and [6|). With the same criterion for the threshold, a 
shallow survey such as PTF will only detect events with 
i^OMopt.off ~ 1, beyond the range of existing on-axis 
SGRB afterglows. In Figure |6] we plot the detectable 
region in the Ej — n phase-space for the maximal LSST 
depth and the LSST normal survey depth (FOM opty oS ^ 
0.01). The allowed phase-space is bounded by the on- 
axis figure of merit (Equation [T]), and we also introduce 
an upper density cut-off of n < 1 cm -3 as an optimistic 
density for a merger in the interstellar medium of a disk 
galaxy. Due to the different dependencies of FOM opit0n 
and FOM opt ,oS on Ej and n, these conditions define a 
triangular region of allowed phase-space for detections of 
off-axis optical afterglows; in §3.21 we perform a similar 
calculation for off-axis radio afterglows to compare the 
relative Ej — n phase-space that is probed by each band. 
Most importantly, we find that the tracks for existing 
SGRB afterglows cross the phase-space region covered 
by optical searches. 

We explore the detectability of on- and off-axis optical 
afterglows more precisely with a Monte Carlo simula- 
tion to determine the fraction of GW events within 200 
Mpc that would be detected by an optical survey with 
a given limiting magnitude and cadence. For targeted 
follow-up searches we use a 1-day cadence with limit- 
ing magnitudes of 22 (PTF), 23.5 (Pan-STARRS), and 
26.5 (LSST); we also include a 1-day and 4-day cadence 
with the standard LSST depth of 24.7 mag. The re- 
sults, summarized in Table [2j show that if the jet energy 
and circumburst density are similar to those required to 
explain the on-axis SGRB data (Figure [3]), then events 
with o bs < 26 j are sufficiently bright to be detected in 
at least 3-5 epochs, given a survey with a depth simi- 
lar to the standard LSST survey (24.7 mag), but with a 
faster cadence of ~ 1 d. Shallower searches are also ca- 
pable of detecting energetic afterglows in a few epochs, 
but this may not be sufficient for a clear identification. 
By contrast, in most cases events viewed at larger angles 
($obs ^ 28 j) are not detectable, even near peak emission 
with LSST. 

The same information is presented graphically in Fig- 
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Fig. 3. — Optic al rest-frame luminosity of existing SG RB afterglows (detections=rect squ ares; upp er limits=fehte triangl es; Bergcrl l2010al : 
Fong ct al. 2011a). Solid lines are afterglow models from van Ecrtcn & MacFadven (2011) (sec also van Ecrtcn ct al. 2010), calculated for 
on-axis observers (0 o b s = ft = 0.2) for a range of jet energies (Ej) and circumburst densities (n). The existing afterglows define an upper 

bound on a figure of merit, FOM op t,on = Ej ^ n^ 2 < 0.1. Also shown are a range of plausible kilonova models (gray shading). The 5a 
limiting magnitudes of various wide-field telescopes are marked by dashed lines; for PTF and Pan-STARRS we assume a maximum of 0.5 
hr per pointing to cover a typical GW error region with a 1-day cadence, while for LSST we show both the normal survey depth and the 
depth for 0.5 hr exposures. 
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Fig. 4. — Same as Figure [3] but for off-axis observers with 
$ob s ~ 26 j . The kilonova emission (gray shading) is isotropic and 
hence remains unchanged for both on- and off-axis observers. The 
range of existing SGRB optical afterglows, covered by the yellow, 
green, and brown lines, indicates that observations with LSST are 
essential. 



ure [7] where we plot contours of detection fraction in 3 
and 5 epochs as a function of depth and cadence. We 
find that in the case of Ej ~ 10 50 erg, the standard 
LSST cadence and depth are sufficient for multiple de- 
tections. However, for lower energies (which may be typ- 
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Fig. 5. — Same as Figuref4]but for off-axis observers with 8 ^, s £S 
48 j. The kilonova emission (gray shading) is isotropic and hence 
remains unchanged for both on- and off-axis observers. The range 
of existing SGRB optical afterglows, covered by the yellow, green, 
and brown lines (well below the limit of the plot), indicates that no 
existing or future telescope will be able to detect optical emission 
at such large off-axis angles. 

ical of most SGRBs), a faster cadence and greater depth 
(~ 26.5 mag) are required for multiple detections. To 
achieve a detection fraction of 50% in 3(5) epochs for the 
case of 6>obs = 20,- requires a depth of at least 23.5(26) 
mag for a 1-day cadence. 
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the tracks for existing SGRB optical afterglow detections (limits) from Figure[3] The existing data suggest that radio detections are highly 
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Since detectable optical emission is limited to off-axis 
angles of < 26 j, we estimate the corresponding fraction 
of GW events with potential optical afterglow detections 
as: 



fa 



Pdct d6 « 6.80f + O 



where 9j <C 1 is the average opening angle, and 

6 cos 



p dDt (9) w 0.152 sinfl (1 



cos 



(3) 



(4) 



is the detection probability of events with inclination an- 
gles between 9 and 9 + d9 (e.g., Schutz 2011; their Equa- 
tion 28). As noted earlier, since this scenario is nearly 
face-on, f opt is a factor ~ 3.4 higher than the detection 

- 2 

fraction »1- cos 29 j « 26 j for isotropic emission. 
_ Equation [3] shows that if the average opening angle is 
6j ~ 0.12, which is the value inferred for G RB 051221A 
(|Burrows et all 120061: iSoderberg et all 120061) . as well as 
the typical opening angle required to reconcile the ob- 
served SGRB rate with the best-bet NS-NS merger rate 
O, then up to / op t ~ 0.1 of GW events will be accom- 
panied by potentially detectable optical afterglows. This 
result is consistent with the rate of a few afterglows per 



year inferred by iCoward et al.l ([201 ID for their assumed 
total ALIGO/Virgo merger rate of ~ 135 yr -1 . On the 
other hand, if 9j is much larger, > , 0.4 (e.g., as found for 
GRB 050724 bv lGrupe et~altl2006l) . then / opt is of order 
unity, but the overall GW event rate may be lower than 
the best-bet ALIGO/Virgo rate. 

Beyond considerations of depth and cadence, a unique 
optical identification of GW events also requires discrimi- 
nation between off-axis afterglows and potential contam- 
inants. We discuss this issue in i}5l 



3.2. Radio Afterglow 

NS-NS/NS-BH mergers may also be accompanied by 
non-thermal radio afterglow emission, which can origi- 
nate either from the ultra-relativistic jet (as in the case 
of the optical aft erglow), or from mo re spherical, sub- 
relativistic ejecta ([Nakar fc Piranll2011t hereafter NP11). 
The latter includes matter ejected dynamically during 
the merger process ("tidal tails"), or in outflows from the 
accretion disk (see Figure [TJ. Adopting standard models 
for synchrotron emission from a relativistic shock, NP11 
estimate that the peak radio brightness for these cases 
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Fig. 7. — Fraction of off-axis optical afterglow events detected in 
3 (solid line) or 5 (dashed line) epochs as a function of the depth 
and cadence of a search. The two models shown - Ej = 10 50 erg, 
n = 10~ 3 cm -3 and Ej = 10 48 erg; n = 1 cm -3 — have the same 
colors as in Figure [3] 
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'„ )P w 40 E jt50 n /8 ^.2 4 difaoo v \ (5) 

where /3q.2 = v c j/0.2c, i/i is the observing frequency in 
GHz; and dh = 200c^ 2 oo Mpc is the luminosity distance, 
again normalized to the ALIGO/Virgo range for NS-NS 
mergers. Equation [5] also assumes characteristic values of 
p = 2.5 for the electron distribution power law index, and 
£e = £b = 0.1 for the fractions of energy density imparted 
to relativistic electrons and magnetic fields, respectively. 
The radio emission peaks at the deceleration time: 
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tde 



^•6 1/3 00.2 
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(6) 



The peak brightness depends sensitively on both the 
properties of the ejecta (E and /3) and on the circum- 
burst density. As we discuss in detail below, the realistic 
detection threshold for a convincing detection with the 
EVLA (even with ~ 30 hr per epoch) is about 0.5 mJy. 
This requirement therefore defines a figure of merit for a 
radio detection of: 

FOM Iad = E jt50 n 7 /s /3 11 / 4 > 0.2. (7) 

With the exception of the velocity parameter, this fig- 
ure of merit is identical to the case of off-axis optical 
afterglows in terms of the dependence on Ej and n. 

For quasi-spherical ejecta, a characteristic mass of 
M j ~ 10~ 2 M Q in tidal tails or disk winds has an en- 
ergy E » Moj^j/2 ~ 10 50 - 10 51 erg for the expected 
range of velocities 10 /? ~ 0.1 — 0.3. This results in at most 
FOAfi a d » 0.4 n^ 8 , requiring no > 1 cm -3 for a detec- 
tion. For mor e typical densities of < 0-1 cm" 3 associate d 
with SGRBs (|Berger et all [20051: [SodeTberg et al.ll200ll . 
the radio emission from quasi-spherical ejecta will be es- 
sentially undetectable unless the energy scale is much 
larger than ~ 10 51 erg (Figure [6]). 

Equations EH3 can also be applied to the case of off- 
axis afterglow emission using (3 « 1 (Na kar fc Fir an 

9 The ejecta m ass and velocity may be higher in some NS- 
BH mergers (e.g., IRosswoil 120051 ) . especially th ose that merge 
on eccentric orbits i n dense stellar clusters (e.g., ILee et all 120101 ; 
IStephens et alJI2011TI . but these are unlikely to represent the typi- 
cal case. 



1201 If) along with values for the jet energy and circum- 
burst density inferred from the optical afterglow data 
(FOM opt ,on < 0.1; Equation P). In Figure M we plot 
the region of E — n phase-space that is accessible to ra- 
dio detections (FOM ra d > 0.2). As can be seen from 
the Figure, none of the existing SGRB optical afterglow 
intersect this region, indicating that radio detections of 
off-axis afterglows are likely to be rare despite the overall 
isotropy of the signal. 

We now address in detail the estimated minimum radio 
brightness necessary for a successful detection. Although 
faint radio emission is in principle detectable with a deep 
integration, a significant challenge is the small field of 
view of sensitive instruments such as the EVLA (w 0.4 
deg 2 at 1 GHz), requiring ~ 100 — 200 pointings to cover 
a typical GW error region of tens of square degrees. Tar- 
geting individual galaxies within the error region does 
not decrease the number of required pointings since there 
are ~ 400 galaxies with L > 0.1 L* within a typical error 
region (to 200 Mpc). Even with only 10 min per point- 
ing, ~ 30 hr per epoch will be required to cover the full 
error region 11 , already a substantial allocation of EVLA 
time. Multiple epochs will be required over a span of 
weeks to years to detect the rise and decline of the radio 
light curve following a GW detection (Equation [6]), for a 
total of about ~ 300 hr of EVLA time (i.e., to search the 
best-bet rate of ~ 40 GW triggers per year will require 
essentially 100% of the EVLA time). Thus, a reasonable 
exposure time per pointing is < 10 min, which at 1 GHz 
corresponds to a 5cr limit 12 of about 0.25 mJy. Since a 
convincing detection will require the brightness to rise to 
about twice the threshold, the minimum detectable peak 
flux is F VtV 0.5 mJy. We note that the threshold may 
be even higher in the compact EVLA configurations (C 
and D) due to substantial source confusion imposed by 
the large synthesized beam size (12 — 44"). 

Observations at a higher frequency of 5 GHz can in 
principle provide better sensitivity (and reduce source 
confusion problems), but in reality will actually require 
even more observing time. This is mainly because the 
field of view at 5 GHz is sufficiently small (0.02 deg 2 ) that 
a more profitable strategy is to target the ~ 400 galaxies 
with L > 0.1 L* within a typical GW error region. Even 
with only 5 min per pointing this will require about 40 hr 
per epoch, with a resulting 5<r limit of 0.1 mJy. A con- 
vincing detection will therefore require F UtP > 0.2 mJy, 
which given a typical spectrum of F v>p oc i/~ - 75 is equiv- 
alent to a limit of > 0.7 mJy at 1 GHz, worse than the 
1 GHz observing strategy, with even more time required 
per epoch. 

Observations with future wide-field radio interferom- 
eters (e.g., ASKAP) will cover a typical GW error re- 
gion with a few pointings, requiring only a few hours per 
epoch. However, these instruments suffer from poorer 
angular resolution compared to what is possible with the 
EVLA (e.g., ASKAP with - 10" resolution). This will 
lead to significant source confusion at the required low 
flux density levels. More critically, radio emission from 
the host galaxy itself will present a challenge; at 200 Mpc 
a star formation rate of only 1 M Q yr _1 corresponds to 

The typical overhead for phase, flux, and bandpass calibration 
with the EVLA is about 25%. 



' https : //science .nrao . edu/f acilities/e via/calibration- and- tools/ex; 
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a 1 G Hz flux density of about 0.6 mJy (|Yun fc Carillil 
12002( 1 . At a resolution of 10" (10 kpc at 200 Mpc) galax- 
ies will generally appear as unresolved point sources, 
and will prevent the detection of significantly fainter co- 
incident radio counterparts. Thus, an instrument like 
ASKAP will cover a GW error region faster than the 
EVLA, but to a similar effective depth limited by source 
confusion. 

A final complication with radio detections is the long 
time delay between a GW trigger and the peak of the 
putative radio signal, which could negate a robust asso- 
ciation. For a sub-relativistic counterpart (0 ~ 0.2) with 
an optimistic density of n ~ 1 cm -3 , a detection requires 
E > 10 51 erg (Equation [7]), and as a result idee ~ 6 yr, 
requiring observations for over a decade. For the rela- 
tivistic case (0 « 1) with n ~ 1 cm~ 3 , the peak time 
corresponding to a detectable signal is idee ~ 0.1 yr. 
The latter case will require a ~ week cadence to ro- 
bustly sample the light curve, corresponding to about 
15 - 20% of the EVLA time (with ~ 30 hr per epoch). 
The absence of a credible detection will require a ~ year 
cadence to search for a non-relativistic counterpart. Of 
course, with a multi-year timescale the probability of 
mis-identification with an unrelated radio transient be- 
comes larger. 

Despite the various difficulties outlined above, a clear 
advantage of radio searches is the lower number of con- 
taminating sources compared to the optical band. As 
discussed in NP11, confusion with AGN radio variabil- 
ity can be reduced by requiring an offset from the cen- 
ter of the host galaxy, although this may be difficult 
with an angular resolution of > 10" (EVLA in its com- 
pact configurations and ASKAP). Similarly, while some 
normal Type Ib/c supernovae have similar radio light 
curves to those expected for NS-NS mergers (since they 
produce ejecta with ~ 0.3), they are generally less 
energetic, with only ~ 10 47 — 10 48 erg coupled to the 
fast ejecta (jBerger et alJl2002L 120031 ) . These events will 
also be accompanied by optical supernova emission on a 
similar timescale, providing an additional source of dis- 
crimination. Finally, relativistic Type Ib/c supernovae 
(with or without an associated GR B) have ~ 10 49 — 10 50 
erg coupled to their fast ejecta (jKulkarni et al.l fl998l : 
ISoderberg et al.l [2010( 1 . but these are also accompanied 
by a bright optical supernovae. 

To conclude, the utility of radio emission as an EM 
counterpart is particularly sensitive to the typical energy 
and circumburst density. In the case of off- axis after- 
glows, detections require a high energy and density that 
exceed those of known SGRB afterglows (Figure [5]) ■ In 
the non-relativistic case, even higher energy and/or den- 
sity are required, such that for an expected upper bound 
of n < 1 cm~ 3 the required energy is E > 10 51 erg. The 
required telescope time for an effective search is hundreds 
of hours (EVLA) , with perhaps only tens of hours using 
future wide-field instruments (e.g., ASKAP). The time 
delays range from months to years, which may compli- 
cate a robust association. The key advantages are the 
spherical geometry at t > idee and the smaller number 
of contaminating sources compared to the optical band. 

4. KILONOVA 

The detectability of SGRBs and their afterglows is sen- 
sitive to uncertainties in the degree of relativistic beam- 



ing and, in the case of afterglows, the properties of the 
circumburst environment. Of course, it is also possi- 
ble that not all NS-NS mergers produce SGRBs. How- 
ever, independent of this association, the mergers are 
expected to be accompanied by isotropic thermal emis- 
sion, powered by the radio active decay of heay y ele- 
ments in t he merger ejecta (iLi fe Paczvri ski 1998; here- 



after LP98; iKulkarnfcboa lRosswodl200 
[20101 : iRoberts et all 120111 : 
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Type la supernovae, which are powered by the decay of 
56 Ni and 56 Co, the ejecta from NS-NS mergers is pri- 
marily neutron-rich (electron fraction Y e <C 0.5) and 
thus produce little nickel. Instead, heavier radioactive 
elements (mass number A > 130) are expected to form 
as neutrons capture onto nuclei (r-process nucleosynthe- 
sis) aft er the ejecta decompresses from nuclear densities 
(e.g.. lLattimer fc Schramm! 11974 lEichler et ail 119891: 
iFreiburghaus et al.lll999f ). "Although the r— process itself 
lasts at most a few seconds, these newly-synthesized ele- 
ments undergo nuclear fission and beta decays on much 
longer timescales. The resulting energy release will power 
bright emission once the ejecta expands sufficiently that 
photons can escape. 

Neutron-rich material is expected to be ejected 
both dynamically d uring the final coalescence (e.g., 
IRosswog et al.l 11999( 1 and b y outflows from the accre- 
tion disk at later times (e.g.. IMetzger et al. 2008a. [20091 : 
iDessart et aTll2009l: [Lee et al.ll2009t Figure [Tj). Depend- 
ing on the properties of the merging binary, expected 
values for the ejecta mass and velocity are in the range 
M ej ~ 10~ 3 - 0.1 M w and ps 1 - .3, respectively 
(e.g. IRosswog et alj[l999t lRosswog||2005D . The resulting 
emission peaks when photons are able to diffuse thr ough 
the ejecta on the expansion timescale (|Arnettlll982l ): the 
low ejected mass thus results in a somewhat dimmer and 
faster evolving light curve than a normal supernova, last- 
i ng days instead of weeks. 

IMetzger et all (|2010f l use a nuclear physics reaction 
network to calculate the radioactive heating of the ejecta 
from NS mergers, and a radiative transfer code to model 
the light curve and color evolution. For typical values of 
M e j = 10~ 2 M and — 0.1, they find that the transient 
peaks at an absolute visual magnitude of My ~ — 15 on 
a timescale of ~ 1 d; because this is approximately one 
thousand times brighter than novae (yet dimmer than 
swper-novae) they dub these ev ents kilonovae 1 3 . 

Although the calculations of IMetzger et al.l (|2010f ) in- 
clude full radiative transfer, they show that the kilo- 
nova light curve is well-approximated using a simple one- 
zone model (LP98), provided that one adopts a value of 
/nuc ~ 3 x 10 -6 for the dimensionless parameter quan- 
tifying the amount of nuclear heating on a timescale 
of ~ 1 d (LP98). Similar resul ts for the rad io active 
hea ting were found rec ently by IRoberts et al.1 (|2011f ) 
and I Gorielv et all (|2011l ). despite somewhat different as- 
sumptions about the geometric structure and thermody- 
namics of the ejecta. 

In Figures [3] and [4] we plot a range of kilonova models 
that span the expected range of ejecta mass and velocity, 
allowing also for realistic theoretical uncertainties in the 
value of / nuc and the opacity of pure r-process ejecta. 

13 The ter ms mini- supernovae (LP98) and macro-novae 
(Kulkarni 2005) are also sometimes applied. 
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Fig. 8. — Fraction of kilonova events detected in 3 (solid line) or 
5 (dashed line) epochs as a function of the depth and cadence of 
a search. We have adopted a representative kilonova model with 
A/ cj = 1(T 2 M©, P = 0.1 c, and / lluc = 3 x 1CT 6 . 

The resulting kilonova emission peaks on a timescale of 
~ 0.5 — 5 d, with an optical luminosity in the range 
~ 10 41 - 10 42 - 5 erg s~\ corresponding to w 19 - 22.5 
mag at the edge of the ALIGO /Virgo volume. Following 
the peak, the kilonova luminosity declines as L v oc t~ a , 
with a ~ 1 — 1.4, due to the declining radioactive power; 
the actual light curve may decline even faster once 7- 
rays or /3-decay leptons freely escape the ejecta without 
depositing their energy. 

An important characteristic of kilonovae are their rel- 
atively unique spectra, which can serve to distinguish 
these events from other astrophysical transients. Overall, 
the kilonova spectrum is predicted to be quasi-thermal 
with T 1=3 10 4 K (although line blanketing in the UV may 
substantially redden the color temperature). Near peak, 
doppler broadening caused by the high ejecta velocity 
will smear out individual spectral features and the overall 
continuum will be smooth. Following the peak, however, 
the photosphere will recede deeper into the ejecta, where 
the velocity is lower. Individual spectral lines from res- 
onant transitions may then become apparent. Since the 
ejecta are composed entirely of exotic heavy nuclei, the 
dominant spectral features may not resemble those of any 
supernova detected to date. Detailed predictions of kilo- 
nova spectra are unfortunately impossible because lab- 
oratory data on the spec tral lines of r-proc ess elements 
are currently sparse (e.g.. lLawler et aLll2009D ; the closest 
known analog to a pure r-process photosp here are ultra 
meta l-poor stars in the Galactic halo (e.g., ISneden et al.l 
l200i . 

To assess the detectability of kilonovae we carry out 
a Monte Carlo simulation for an optical surveys with 
a range of limiting magnitudes and cadences (Table [3]). 
We find that at the depth and cadence of the normal 
LSST survey (r ss 24.7 mag, At = 4 d), essentially no 
kilonovae will be detected in 5 epochs, unless M e j ~ 0.1 
M Q ; about 3/4 of all events will be detected in 3 epochs 
with the normal LSST survey if M oj ~ 10" 2 M Q . To 
detect events with M e j ~ 10 -2 M Q in 5 epochs requires 
a 1-day cadence, preferably with telescopes capable of 
reaching > 23 mag (e.g., Pan-STARRS, LSST). Finally, 
for M j ~ 10 -3 M©, no existing or planned telescope will 
provide 5 detections, but LSST with a 1-day cadence is 
likely to provide 3 — 4 detections. These results are also 



summarized in Figure [8j where we plot contours for the 
fraction of kilonovae detected in 3 and 5 epochs as a func- 
tion of limiting magnitude and cadence, assuming typical 
values of M oj = 10~ 2 M©, B = 0.1, and / nuc = 3 x 10" 6 . 
The plot demonstrates that to achieve 50% completeness 
in 3(5) epochs given a cadence of ~ 1 d requires a limit- 
ing magnitude of > 21(22.5) mag. We discuss potential 
contamination from other optical transients in fj5] 

5. SUMMARY: WHAT IS THE MOST PROMISING EM 
COUNTERPART? 

We now bring together our conclusions from the previ- 
ous sections to address the question of the most promis- 
ing EM counterpart. A summary of the expected de- 
tection fractions, dependence on density, and the Cardi- 
nal Virtues satisfied for each EM counterpart is provided 
in Table Q] We first discuss the case in which all NS- 
NS mergers are accompanied by SGRBs, and hence by 
on-axis 7-ray/afterglow emission or by off- axis afterglow 
emission; we then turn to a discussion of the kilonova- 
dominated case. 

5.1. Gamma-Rays 

Short GRBs are easily detectable within the 
ALIGO/Virgo volume with current 7-ray satellites in 
cases when 9 h s < 0j\ they therefore satisfy Virtue #1. 
Although this configuration applies to only a small frac- 
tion of all mergers (and therefore violates Virtue #2), 
the SGRB rate within the ALIGO/Virgo volume (en- 
hanced by a factor of 3.4 for face-on mergers) is suffi- 
ciently high that ~ 1 coincident event should occur per 
year (Figure (2). SGRBs thus represent an ideal counter- 
part to confirm the cosmic origin of at least some GW 
events, and to test whether SGRBs in fact accompany 
NS-NS/NS-BH mergers. Such an association is critical 
since it will help to justify the expensive search for or- 
phan afterglows in the optical and radio bands. SGRBs 
also suffer from little contamination and therefore sat- 
isfy Virtue #3. It is therefore critical that a sensitive 
7-ray satellite be in operation during the ALIGO/Virgo 
era. Fine positional accuracy for an SGRB detection 
(e.g., Swift) is less critical than all-sky coverage (e.g., 
Fermi/ 'GBM) since the temporal association alone within 
the large error region of a GW source would suffice to 
determine an association. A ~ arcsecond position (sat- 
isfying Virtue #4) could then be achieved from the ex- 
pected on-axis optical afterglow or a kilonova, which will 
be brighter than ~ 22 mag and hence easily detectable 
with wide-field telescopes. Thus, SGRBs satisfy 3 out of 
the 4 virtues for a promising EM counterpart. 

5.2. Off-Axis Optical and Radio Afterglows 

In the absence of 7-ray emission, orphan afterglow 
emission (both optical and radio) is the most promis- 
ing counterpart if the typical jet energy and circum- 
burst density lie near the upper end estimated from cur- 
rent SGRB observations: Ej^ n^ 2 ~ 0.1. For off-axis 
optical afterglows the detectability limit with LSST is 

FOM optioS = E^ 50 n /8 > 0.002, as long as 9 ohs < 26 j. 
Thus, optical afterglows satisfy Virtue #1, but violate 
Virtue #2 since a fraction of at most ~ 79? ~ 0.1 
would be detectable. For radio afterglows the detectabil- 
ity limit with EVLA or a future instrument like ASKAP 
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TABLE 1 

Comparison of Electromagnetic Counterparts 



Counterpart 


Detection Efficiency 


Depends on 


Virtues 


Follow-up 




Fraction 


Density? 


Satisfied 


Instruments 


Short GRB ({gj 


~ 3.4 x f bn x FOV 7 a 


no 


1,3 


Fermi/GBM 


Orphan optical afterglow (i|3.1l) 


~ 78? x J r o P t(E i ,n) b < 0.1 


yes 


1,3(?),4 


Pan- STARRS, LSST 


Orphan radio afterglow f i|3,2H 


- 1 x •7 r r ad(£j,«) c 


yes 


1,3,4(?) 


EVLA, ASKAP 


Non-relativistic radio f£|3.2ll 


?; only if E > 10 51 erg 


yes 


1(?),3(?),4(?) 


EVLA, ASKAP 


Kilonova (0) 


~ l d 


no 


1,2,3(?),4 


PTF, Pan-STARRS, LSST 



Note. — a Field of view of gamma-ray telescope as a fraction of 47r steradian; fraction of mergers accompanied by a sufficiently 
energetic jet and dense circumburst medium for an optical detection (related to FOM opt o ff = Ej ^n 7 ^ & > 0.002; Equation [2t: 
c fraction of mergers accompanied by a sufficiently energetic jet and dense circumburst medium for a radio detection (related to 
FOM ra d = Ej ^Q n^/ 8 > 0.2; Equation[7J; d assuming that a telescope similar to LSST covers the GW sky error region with a 
cadence of ~ 1 day. 



is FOATrad = -Ej-,50 nj/ 8 > 0.2 independent of viewing 
angle (i.e., they satisfy Virtue #2). However, as a result 
of the limited range of E — n phase-space probed by ra- 
dio observations, they violate Virtue #1. Indeed, using 
the tracks for existing SGRB optical afterglows in the 
E — n phase-space, we find that none cross the portion 
accessible to radio searches (Figure [6} . It is therefore 
possible that despite the relative isotropy of the radio 
emission, existing and planned instruments are simply 
not sensitive enough to detect the emission for typical 
SGRB parameters. On the other hand, about half of all 
existing optical afterglows will be detectable to the depth 
of LSST (Figure O, but only with a viewing angle up to 
~ 28 j, indicating an expected optical detection fraction 
of<5%. 

In terms of Virtue #3, contamination in the radio 
band is less severe than in the optical band. In the 
optical, we expect contamination mainly from back- 
ground supernovae (Type la, and Type IIP shock break- 
out), while contamination from AGN variability can be 
avoided based on its coincidence with a galaxy nucleus. 
We evaluate the off-axis light curves in comparison to 
supernova light curves in Figure H] For the high density 
case (n ~ 1 cm -3 ), the off-axis light curve peaks on a 
timescale of ~ 1 d, and is indeed similar to the kilonova 
light curves. The rapid rise and decline (declining by 
~ 3 mag in ~ 5 d) can easily distinguish this case from 
Type la supernovae, which rise and decline by compa- 
rable amounts on timescales of tens of days (Figure^. 
Shock breakout emission from red supergiants also leads 
to rapid rise and decline (~ 1 mag on a timescale of ~ 1 
d), but is subsequently followed by a long and bright 
plateau phase (a Type IIP supernova) that can easily 
distinguish these cases. Thus, with a sufficiently rapid 
cadence (1-day) and a depth similar to the LSST nor- 
mal survey (or better yet ~ 0.5 hr pointings with sa 26.5 
mag), off-axis afterglows in a dense medium can be sep- 
arated from background contaminating supernovae. 

The case of an off-axis afterglow in a low density 
medium (n ~ 10~ 3 cm~ 3 ) is somewhat more compli- 
cated. The off-axis light curve peaks on a timescale of 
~ 20 d, followed by a decline of about 1 mag in the sub- 
sequent 3 months (Figure [9|) . With a peak brightness of 
~ 24.5 mag, a convincing detection requires a depth be- 
yond the normal LSST survey mode. However, a 1-day 
cadence is not essential, and the depth can be achieved by 
stacking multiple images ona~ week timescale. Given 
the slower evolution of the light curve, it is more similar 



to supernova light curves than the n ~ 1 cm -3 case. The 
predicted rate of decline is slower than a Type la super- 
nova post maximum (i.e., « 2.5 — 3 mag in ~ 100 d). 
It is, however, faster than a typical Type IIP supernova 
light curve, which exhibits a plateau for ~ 100 d. 

Color evolution can in principle also be used to dis- 
tinguish off-axis afterglows from other transients. Since 
the afterglow is synchrotron emission (and the optical 
waveband is generally above the characteristic frequency, 
v m )i it has a power-law spectrum with a fixed slope 
F v oc j/( 1_ p)/ 2 and hence a constant red color g — r~ 0.2 
mag for p = 2.5. By comparison, Figure [9] shows that 
the color of a shock break-out and rising IIP SNe in- 
creases by a magnitude from blue to red in just a few 
days. Although the colors of a rising SN la are similar to 
the afterglow emission, events observed near their peak 
(the case in which a background Type la supernova light 
curve could be mistaken for a low density afterglow) are 
much redder. 

Finally, both radio and optical counterparts will satisfy 
Virtue #4, although at low frequency and low signal-to- 
noise ratio, EVLA/ ASKAP positions will typically be 
> few arcsec, as opposed to sub-arcsecond in the optical 
band. At a typical distance of 200 Mpc this should not be 
an impediment for a host galaxy association (1" ss 0.8 
kpc), but it will not allow a robust study of the sub- 
galactic environment, and he nce an association with spe- 
cific stellar populations (c.f., iFong et "all I2010I ). It may 
also impede the rejection of AGN. 

We therefore conclude that optical and radio afterglows 
do not satisfy all of the required cardinal Virtues for an 
EM counterpart. The fraction of detectable off- axis opti- 
cal afterglows is ~ 0.1, and possibly even lower depending 
on the range of energy and circumburst density for typ- 
ical NS-NS/NS-BH mergers. The fraction of detectable 
radio afterglows may be close to zero due to the limited 
range of E — n phase-space accessible with existing and 
planned radio telescopes. 

5.3. Kilonova 

If the majority of NS-NS/NS-BH mergers occur in low 
density environments (n < 10~ 3 cm~ 3 ) or produce low 
energy jets (E < 10 49 erg), then optical afterglows are 
no longer effective counterparts. This is also true if most 
NS-NS mergers are not accompanied by SGRBs. In these 
cases, kilonovae provide an isotropic source of emission 
that does not depend on the external environment. Be- 
cause the emission is thermal and requires only a small 
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Fig. 9. — Comparison of off-axis afterglow and kilonova light curves with contaminating optical transients, including r-band light curve 
(top-left), V-band light curve (top-right), V — r color curves (bottom- left), and r-band decline rate (bottom-right). Off-axis light curves are 
shown for the parameters in Figures I3H51 We use a typical model for the kilonova light curve, Af ej = 10" 2 Mq and /3 = 0.1, for two different 
assum ptions about t he opacity of the ejecta: the LP98 model assuming black body emission (purple), and a model assuming pure Fe opacity 
from lMetzger et al.l IpOlOT ) (blue ). Shown for comp arison are background Type la supernova light curves 12 d before maximum (black) an d 
1 d after maximum (cyan) from Wang ct al. (2009J), as well as supernova shock break-out from a red supergiant (red) (Nakar & Sari 2010), 
followed by a standard Type IIP supernova plateau (dashed red). 



quantity of neutron-rich ejecta (as is likely to accompany 
most mergers), the predicted signal is also relatively ro- 
bust, and at a peak optical brightness of ~ 19 — 22 mag 
is detectable with a facility such as LSST. Thus, kilono- 
vae satisfy Virtues #1, #2, and #4, as long as a rapid 
and deep search is carried out (Figure [8] and Table [3|). 
Therefore, a key question is whether these events can be 
easily distinguished from contaminating sources. 

We present kilonova light curves in comparison to back- 
ground supernovae in Figure [5] We show both a black 
bod y model (LP98) and a model assuming pure Fe opac- 
ity (|Metzger et al.ll2010f ) to span the plausible range in 
the true light curve and color evolution, the latter of 
which remains especially uncertain due to the lack of ex- 
perimental data on the opacity of pure r-process ejecta. 
As in the case of off-axis afterglow emission in a dense 
medium, the kilonova light curve evolution is much more 
rapid than for supernovae: the rise time is ~ 1 d, followed 
by a decline of about 3 mag in ~ 5 — 8 d. This behavior 
places stringent constraints on an effective search (1-day 
cadence and a depth of > 24 mag), but it allows for a 
clean separation from contaminating supernovae. Thus, 
it appears that kilonovae can satisfy all four Cardinal 
Virtues. 



5.4. Quality of Information: A Fifth Virtue? 

Even if all types of EM counterparts discussed in this 
paper will eventually be detected in conjunction with 
GW triggers, each provides distinct information about 
the merger. The detection of an SGRB in temporal co- 
incidence with a GW trigger will establish a firm connec- 
tion with NS-NS/NS-BH mergers. Since these events are 
also expected to be face-on mergers, such a detection will 
establish the orientation of the binary, thereby allowing 
for more accurate extraction of additional binary param- 
eters from the GW signal, such as the masses and spins of 
its members and, potentially, information about the high 
density equation of state. While a 7-ray detection is not 
itself likely to substantially improve the positional accu- 
racy, an on-axis optical afterglow, or a kilonova, should 
be detectable and will provide a host galaxy association 
and redshift. 

The detection of an orphan optical afterglow (on- or 
off-axis) will provide much of the same information. 
Namely, it will establish a connection with SGRBs and 
will also possess a nearly face-on orientation (since or- 
phan optical afterglows are only detectable to < 29 j). 
The brightness of the optical emission will also provide 
information on the combination of energy and circum- 
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burst density. A radio detection will provide no informa- 
tion on binary orientation (since essentially all off-axis 
angles are detectable). The peak time and brightness 
will provide information on the combination of energy, 
density, and ejecta velocity. As a result of the velocity 
degeneracy a radio detection will not necessarily estab- 
lish an association with SGRBs (i.e., the specific case of 
j3 rs 1). Finally, kilonova detections provide a unique 
probe of the inner workings of the merger, since their 
light curves depend on the mass, velocity, and geometry 
of the ejecta, while their opacity and spectral features 
probe the ejecta composition (spectroscopy will require 
real-time identification near peak). The discovery of a 
kilonova event will also represent the first in-situ obser- 
vation of freshly produced r-process material, the origin 
of which remains perhaps the biggest mystery in nuclear 
astrophysics. However, such detections will not help to 
establish a connection with SGRBs. 

To summarize, the potential connection of NS-NS 
mergers with SGRBs provides a useful EM counterpart 
for both on- and off-axis emission. However, in essen- 
tially all possible scenarios, only a small fraction of GW 
events (< 10%) will be followed by a detectable SGRB 
or optical/radio afterglow. On the other hand, isotropic 
kilonovae will likely provide a larger detection fraction, 
as long as the typical ejected mass is > 1CP 3 M and 
deep observation with 1-day cadence are carried out. 

6. FOLLOW-UP STRATEGY RECOMMENDATIONS 

Taking advantage of any of the potential EM counter- 
parts discussed in the previous sections requires a careful 
observing strategy. In this section we make specific rec- 
ommendations for searches at 7-ray, optical, and radio 
wavelengths to enhance the detection probability, given 
a reasonable allocation of resources. 

6.1. Gamma-Rays 

In the case of 7-ray detections of an associated SGRB, 
the limiting factor is the small fraction of on-axis events 
within the ALIGO/Virgo detection volume. The on- 
axis orientation provides a boost to the detection vol- 
ume, leading to a detectable rate perhaps as high as ~ 1 
event per year. As a result of this low rate, all-sky cov- 
erage in 7-rays is more critical than the ability to sub- 
stantially refine the GW positions. This will still allow 
for a robust association due to the temporal coincidence, 
and the positional refinement can then be achieved from 
the bright on-axis optical afterglow. In addition, since a 
temporal coincidence with an SGRB allows for a lower 
threshold GW detection, a strategy of searching the GW 
data stream based on 7-ray triggers may actually lead 
to a larger number of detections than the opposite ap- 
proach (i.e., it will boost the accessible volume by more 
than a factor of 3.4 times for face-on orientation). Thus, 
we strongly recommend an operational 7-ray satellite in 
the ALIGO/Virgo era, with capabilities similar to the 
Fermi/GBM. 

6.2. Optical 

The search for orphan optical afterglows requires wide- 
field telescopes capable of achieving a depth of at least 
~ 23 mag, and perhaps ~ 26.5 mag for detections of 
typical events. The maximum achievable depth is de- 
termined by the need to cover tens of square degrees 



with a 1-day cadence. Since ~ 10 pointings are required 
with facilities such as PTF, Pan-STARRS, and LSST, 
the maximum time per pointing is about 0.5 hr, leading 
to depths of ~ 22 mag for PTF, ~ 23.5 mag for Pan- 
STARRS, and - 26.5 mag for LSST. For a localization 
region of a few square degrees (i.e., a single pointing), the 
achievable limiting magnitudes are about 1.2 mag deeper. 
The results of Monte Carlo simulations of the detection 
fractions in 3 and 5 epochs for on- and off-axis afterglows 
given a limiting magnitude and cadence are summarized 
in Table [2] and Figure [7] We find that only in the case 
of E ~ 10 50 erg the standard LSST cadence and depth 
are sufficient for multiple detections. However, for lower 
energies (which may be typical of most SGRBs), a faster 
cadence and greater depth (~ 26.5 mag) are required for 
multiple detections. To achieve a detection fraction of 
50% in 3(5) epochs for the case of Ej ~ 10 48 erg and 
o bs = 26*j requires a depth of at least 23.5(26) mag for a 
1-day cadence. Thus, we conclude that standard LSST 
depth and cadence are non-ideal for detections of off-axis 
afterglows. 

A key issue discussed in previous papers is that at the 
typical distance limit of ~ 200 Mpc for ALIGO/Virgo de- 
tections, one could expedite the search for EM counter- 
parts by focusing on galaxies within this volume. How- 
ever, the number of galaxies within a typical GW error 
region in the typical luminosity range of SGRB hosts 
(L > 0.1 L*) is 14 ~ 400; the number of L* galaxies is 
about 50. Thus, the number of galaxies is much larger 
than the required number of pointings for wide-field tele- 
scopes (~ 10), and therefore a focus on nearby galaxies 
has no effect on the required cadence and depth. Con- 
versely, the number of galaxies is too large for an efficient 
search with a large-aperture but small field-of-view tele- 
scope (e.g., Keck, Gemini) since several hundred point- 
ings will be required within a single night (leading to < 1 
min per pointing) . Even a search of only L > L* galaxies 
(which will inevitably miss a substantial fraction of the 
counterparts in sub-L* galaxies) will limit the observa- 
tions to ~ 5 min per galaxy, and will require the full use 
of an 8-m class telescope for several nights. 

Thus, our key recommendation is follow-up with wide- 
field optical telescopes capable of reaching a depth of > 
23 mag in 0.5 hr, using a 1-day cadence. Effectively, this 
means that LSST should execute a non-standard cadence 
to follow up GW triggers (a "sub-survey" mode). Ideally, 
such dedicated follow-up observations will also involve 
longer exposure times than the normal survey mode (up 
to ~ 0.5 hr per pointing), but even without a change to 
the standard exposure time, repeated visits on a nightly 
basis will provide the most efficient search strategy for 
optical counterpart searches. The same strategy is key 
for detections of the fast-evolving kilonovae (Figure[8]and 
Table [3|). If no convincing counterpart is detected within 
a few days, a search for delayed off-axis emission (due to 
low density) can employ the normal LSST cadence since 
the typical timescale is tens of days (Figure ■ 

We note that efforts to perform GW-triggered optical 
follow-up have already begun during the recent LIGO 

14 We use the SDSS luminosity function with M* ~ —21.2 mag 
and (f>* M 5 X 10 -3 Mpc -3 (Blanton ct al. 2003). Integration down 
to 0.1 L* therefore gives about 8 galaxies per square degree within 
a distance of 200 Mpc. 
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science run by the LOOC UP (Locating and Observ- 
ing Optical Counterparts to Unmodeled Pulses) project 
dKanner et alJl200l lAbbott et alJl200l . LOOC UP re- 
constructs the sky position of candidate GW signals to 
make prompt optical follow-up ob servations usin g wide- 
field sub-meter class telescopes. ICannon et al.l ()2011[ ) 
discuss techniques to reduce the latency between GW 
detection and follow-up to minutes or even seconds, in 
which case counterpart searches could in principle be- 
gin simultaneous with the final coalescence. The meter 
class telescopes employed in LOOC UP were sufficiently 
sensitive to detect off-axis afterglows or kilonova within 
the LIGO volume, but will be clearly insufficient for the 
ALIGO/Virgo volume. 

6.3. Radio 

The search for radio counterparts is complicated by the 
highly uncertain peak time of the light curve (ranging 
from months to years), as well as by the small field-of- 
view of the EVLA. For an off-axis afterglow in a dense 
medium the expected peak time is ~ 0.1 yr, while for a 
non-relativistic spherical counterpart the expected peak 
occurs on a timescale of a few years. As a result, a robust 
search has to cover a broad range of timescales, with an 
initial rapid cadence of ~ 1 week, followed by a transition 
to monthly and then yearly observations. In total we 
estimate that at least 10 — 15 epochs will be required, 
spread logarithmically over a decade. As discussed in 
£13.21 a search with the EVLA will require about 30 hr 
per epoch to cover a typical GW error region. If we 
assume that about half of these epochs will take place 
within the first year after the trigger, the required time 
to follow up the best-bet rate of 40 GW triggers per year 
is about 100% of the EVLA time. 

The peak flux density required for a convincing detec- 
tion (a factor of 2 rise above the achievable 5er = 0.25 
mJy threshold) is about 0.5 mJy As shown in Figure [HI 
none of the existing SGRB optical afterglows will lead 
to such a bright signal. The detection of events with 
F u ,p ~ 0.1 mJy is complicated not only by the excessive 
amount of required telescope time (hundreds of hours per 
epoch!), but also by source confusion at low frequency, 
and the fact that a host galaxy with SFR ~ 1 M Q yr -1 
has a 1 GHz flux density of ~ 0.6 mJy. These factors 
will limit the searches even if the error regions are only 
a few square degrees. 

In terms of the search strategy, at 1 GHz the best ap- 
proach is to tile the full error region with ~ 100 — 200 
pointings, while at higher frequencies it is more profitable 
to target the ~ 400 galaxies with L > 0.1 L* individually 
(although this will require even more observing time). 
Observations with a future facility such as ASKAP will 
reduce the time requirement to a few hours per epoch 
by reducing the number of pointings. However, the sen- 
sitivity of the search is unlikely to improve since source 
confusion becomes a dominant obstacle. 

We finally note that for the case of radio emission from 
non-relativistic ejecta, the required energy and density 
are E > 10 51 erg and n > 0.1 cm~ 3 (Figure |6|). It is 
unclear if the required energy scale can be produced in a 
typical NS-NS merger, but even if it does, the resulting 
decade-long delay between the GW trigger and peak of 
the putative radio emission will require many observing 
epochs, and will furthermore impede a convincing asso- 



ciation. 

Thus, our recommendation in the case of radio searches 
is to limit the search to a timescale of a few months, ap- 
propriate for the case of an off- axis afterglow [fi ~ 1). 
In this case the peak brightness is also more likely to be 
detectable, since F„ tP cx /3 11 / 4 , and the relatively modest 
time delay (idee °c /3~ 5 ^ 3 ) will reduce the potential for 
contamination. The required observing time with the 
EVLA will be about ~ 200 hr for several epochs loga- 
rithmically spaced in time over several months. Future 
searches with facilities like ASKAP should also focus on 
the same timescale since they are limited to the same 
depth as EVLA. 

7. CONCLUSIONS 

With the era of gravitational wave astronomy fast ap- 
proaching we investigated and critically assessed a range 
of potential electromagnetic counterparts for NS-NS/NS- 
BH mergers, and their detectability with existing and up- 
coming telescopes. We used the rates of (on-axis) SGRBs 
to predict the detection rate with an all-sky 7-ray mon- 
itor, and existing information about SGRB afterglows 
to predict the appearance and brightness of off-axis op- 
tical and radio emission. Finally, we assessed the light 
curves and detectability of kilonovae. Overall, we found 
that none of the potential EM counterparts is guaran- 
teed to satisfy all four Cardinal Virtues — detectability, 
high fraction, identifiability, and positional accuracy — 
but that critical insight into the merger physics can be 
gained from any of the counterparts for at least some 
events. In particular, we found that: 

• Gamma-ray observations are critical for establish- 
ing a firm connection between SGRBs and NS- 
NS/NS-BH mergers. Such detections are likely to 
be limited to a rate of < 1 yr -1 , but the face-on 
configuration will allow for better measurements of 
the binary parameters. In addition, 7-ray triggered 
GW searches may enhance the probability of joint 
GW/EM detections. 

• The number of expected galaxies with L > 0.1 L* 
(typical of SGRB hosts; iBergerl |2009| ) in a typi- 
cal GW error region is ~ 400, making targeted 
searches of galaxies in the optical and/or radio in- 
efficient. In both cases, complete coverage of the 
error region is less time consuming. 

• On-axis optical emission typical of existing SGRB 
afterglows is easily detectable with existing and 
planned wide-field telescopes at < 200 Mpc. 

• Off-axis optical afterglow emission is only de- 
tectable to o bs ~ 20j, and is hence limited to 
< 10% of all mergers. Within this range, LSST 
observations are required to detect events similar 
to existing SGRB afterglows, but with a special- 
ized depth/cadence of ~ 26.5 mag (achievable in 
0.5 hr) and 1 d. Observations with the normal 
LSST survey mode are likely to miss most coun- 
terparts. With our proposed LSST depth/cadence 
contamination from other sources can be avoided 
based on the rapid rise and decline time for high 
density cases, and based on the decline rate and 
lack of color evolution for low density cases. 



Electromagnetic Counterparts of Neutron Star Mergers 



15 



• Off-axis radio afterglow emission can in principle 
be detected at all observer angles, but existing 
and planned telescopes are limited to flux levels 
that exclude the detection of existing SGRB after- 
glows. The long delay relative to the GW trig- 
ger (up to many years) provides an additional ob- 
stacle to robust association. Still, searches with 
a weekly-monthly cadence and per-epoch exposure 
times of ~ 30 hr (EVLA) or - few hr (ASKAP), 
may lead to detections of rare energetic events in a 
dense medium. Contamination from other sources 
is less severe than in the optical band, but the posi- 
tional accuracy of a detection will be poorer than in 
the optical, potentially preventing studies of sub- 
galactic environments. 

• Isotropic optical emission powered by the radioac- 
tive decay of r-process elements in the merger 
ejecta (kilonova) is expected to reach a peak bright- 
ness of ~ 19 — 22 mag at At ~ 1 d, with a subse- 
quent decline by several magnitudes in a few days. 
The brightness is independent of the ambient den- 
sity. The detection of kilonovae therefore requires 
rapid cadence (~ 1 d) to at least the normal LSST 
survey depth (and preferably a maximal depth of 
~ 26.5 mag). The rapid rise and decline will re- 
duce the contamination from other optical tran- 
sients (e.g., supernovae). With the wherewithal to 
carry out such a search with LSST, kilonovae can 
indeed serve as the most promising counterpart of 
compact object binary mergers. 

• Our key recommendations for maximizing the de- 
tection probability of EM counterparts are: An all- 
sky 7-ray satellite similar to Fermi/GBM; a spe- 
cialized LSST "sub-survey" mode with a 1-day ca- 
dence and a depth of ~ 26.5 mag; radio follow-up 
with a weekly cadence to a depth of 0.25 mJy using 
EVLA/ ASKAP, limited to < few months after the 
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trigger. 



Since the timescales of the various potential EM coun- 
terparts are spread from seconds (7-rays) to days (opti- 
cal) to months (radio) , a staggered approach will clearly 
inform a joint observational strategy. For example, the 
detection of 7-ray emission should trigger an immediate 
high-cadence search for on- axis afterglow/kilonova emis- 
sion in both the optical and radio; such a search will still 
require wide-field imaging. Similarly, the detection of 
off- axis afterglow or kilonova optical candidate (s) should 
trigger targeted radio observations (allowing for much 
deeper observations relative to a complete search of the 
GW error circle). 

We finally note that our framework for evaluating po- 
tential EM counterparts can be revised as the actual posi- 
tional capabilities of ALIGO/ Virgo and the merger rate 
become clear. However, unless the GW sensitivity or 
rate have been substantially over-estimated, we conclude 
that a concerted follow-up effort will determine whether 
NS-NS/NS-BH mergers are associated with SGRBs, and 
will provide critical insight into the physics of compact 
objects and the merger process. The fundamental im- 
portance of these results justifies the proposed expensive 
electromagnetic observational strategy. 
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Detection Efficiency of Orphan Afterglows within 200 Mpc 
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magnitude; (ii) observing cadence; (iii) jet energy; (iv) density; (v) 
viewing angle; (vi) fraction of events detected in one epoch; (vii) frac- 
tion of events detected in three epochs; and (viii) fraction of events 
detected in five epochs. 
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TABLE 3 

Detection Efficiency of Kilonovae within 200 Mpc 
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Note. — The columns are (left to right): (i) 5cr r-band limiting mag- 
nitude; (ii) observing cade nce; (in) ejecta mass; (iv) ejecta velocity; (v) 
nuclear heating parameter (j Li fc P aczvhski 1998); (vi) fraction of events 
detected in one epoch; (vii) fraction of events detected in three epochs; 
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